What do gas-rich galaxies actually tell us about modified Newtonian dynamics? 
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It has recently been claimed that measurements of the baryonic Tully-Fisher relation (BTFR), 
a power-law relationship between the observed baryonic masses and outer rotation velocities of 
galaxies, support the predictions of modified Newtonian dynamics for the slope and scatter in the 
relation, while challenging the cold dark matter (CDM) paradigm. We investigate these claims, and 
find that: 1) the scatter in the data used to determine the BTFR is in conflict with observational 
uncertainties on the data; 2) these data do not make strong distinctions regarding the best-fit BTFR 
parameters; 3) the literature contains a wide variety of measurements of the BTFR, many of which 
are discrepant with the recent results; and 4) the claimed CDM "prediction" for the BTFR is a 
gross oversimplification of the complex galaxy-scale physics involved. We conclude that the BTFR 
is currently untrustworthy as a test of CDM. 

PACS numbers: 95.35.+d, 04.50.Kd, 98.56.Wm 



Introduction. — Despite the plentiful evidence for the 
existence of dark matter in the Universe (see Ref. [1] 
for a recent review), the modified Netwonian dynamics 
(MOND) hypothesis — that Newtonian gravity departs 
from its expected behaviour below a certain acceleration 
scale, thus potentially eliminating the need for any non- 
luminous, non-baryonic matter [2] — has persisted since 
its proposal almost three decades ago. It has recently 
been claimed by McGaugh in Ref. [3] (hereafter MG11) 
and in a follow-up paper [4] that measurements of the 
baryonic Tully-Fishcr relation (BTFR) [5, 6], a scaling 
relation between the baryonic masses and outer rotation 
velocities of galaxies, match precisely with a prediction 
of MOND involving zero free parameters, with scatter 
about this relation attributable solely to observational 
uncertainties. MG11 also claims that the observed BTFR 
deviates significantly from that predicted by the stan- 
dard cold dark matter (CDM) scenario. Investigations 
into each of these claims reveal flaws that substantially 
weaken the conclusions of MG11, and highlight the diffi- 
culty of using such an apparently simple relationship as 
the BTFR as a clean test of new physics. 

The data for the 47 galaxies used in MG11 are collected 
from three sources [7-9], subject to the criteria that the 
mass of molecular gas in each galaxy exceeds its stellar 
mass, and that each galaxy has a resolved rotation curve 
that asymptotes to a constant velocity at large radius. 
We have assembled our own sample of galaxies (hence- 
forth FS11) from the same three sources, subject to the 
same two criteria and re-calculating derived quantities 
where needed. 

We find 58 galaxies that are suitable for use. In 
brief, we include 11 galaxies from Ref. [7] and one from 
Ref. [8] that we find to meet the necessary criteria but 
that are excluded from MG11, and exclude one galaxy 
used by MG11 that we find to have star-dominated 
mass. Many of the extra galaxies have low inclina- 
tions (i < 45°), which have been seen to decrease the 



slope of the BTFR [7] , but we see no clear quality issues 
that should cause these data should be excluded (espe- 
cially since MG11 utilizes galaxies with inclinations as 
low as 29°). 

However, we should note that the criteria of MG11 
likely impart strong selection effects and biasing in the 
determination of the BTFR. For example, as discussed 
by Ref. [10], the observability of a flat region in a rota- 
tion curve favours galaxies for which the distribution of 
neutral hydrogen is more extended compared to the ra- 
dial scale of the outer halo — indeed, many otherwise well- 
behaved galaxies do not possess fiat rotation curves [11]. 
In addition, galaxies undergoing significant interaction 
with others nearby will tend to be disfavoured. Fur- 
ther criteria based on data quality are sometimes used in 
BTFR studies [4] , but imposing such criteria can tend to 
skew the sample towards galaxies with particular proper- 
ties, when in fact the goal should be to use a broad range 
of data to avoid selection affects biasing the results. 

Scatter in galaxy data. — The strategy of MOND is 
most simply characterized by its suggestion that the ac- 
celeration on of a test particle in Newtonian gravity is 
actually related to the observed acceleration a by the re- 
lation 



a N = /x(a/a )a, 



(1) 



where ao is a new physical constant with dimensions of 
acceleration and magnitude roughly 10 -10 ms -2 , and 
n(x) is a smooth function that asymptotes to unity for 
x ^> 1 and to x for i <C 1. The modern perspective is 
that the relation (1) might arise from some fully covariant 
theory of non-Einsteinian gravity [12-14], although such 
theories continue to face difficulties of their own [15-18]. 

Upon taking the acceleration of a test particle within 
a galaxy to satisfy a <C do, and assuming that galaxies 
contain only baryonic matter, simple algebra yields the 
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FIG. 1: Histograms (green bars) of values of a — V* /GM^ (in 
units of ms -2 ) for MG11 and FS11 data sets, along with fitted 
Gaussians (black, solid curves) and Gaussians with the most 
likely mean and variance for data sets with the same statistical 
properties as the measured data (blue, dashed curves). For 
MG11, the Monte Carlo-generated distribution is much wider 
than the fitted distribution, indicating a discrepancy in the 
level of scatter in the data as compared to the stated error 
bars. This discrepancy is less severe for FS11. 



MOND prediction for the BTFR: 



Ga ' 



(2) 



where Mb is the (baryonic) mass of a galaxy and Vf is 
the rotation speed asymptotically approached at large 
radius. For disk galaxies, clq should be multiplied by 
an extra factor \ « 0.8 to account for the difference 
in rotation speeds between disk-like and spherical mass 
distributions [4]. A good indicator of the scatter of a 
sample of galaxies about this relation is the distribution 
of a values derived from the measured values of Mb and 
Vf for the sample. 

Therefore, analogously to MGll's figure 3, we calcu- 
late a = V f 4 /GM\y for each of the 47 galaxies used by 
MG11. In Fig. 1, we display a histogram of these values, 
as well as a Gaussian fit to this histogram. The fitted 
Gaussian has width a = 0.28 dex, similar to the width 
of the Gaussian from MG11 (which was not fitted, but 
rather estimated by eye), quoted as 0.24 dex. 

To put the scatter in a values in perspective, we gen- 
erate a large number (10 7 ) of simulated datasets of 47 
galaxies [i.e. (Vf,Mb) pairs], with the velocity and bary- 
onic mass of each galaxy drawn from a normal distri- 
bution determined by the error bars on the correspond- 
ing galaxy from the MG11 sample. We calculate the 
mean and variance of the a values for each dataset, and 
from this determine the most likely mean and variance 
of a galaxy sample with the same properties as that of 
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FIG. 2: 68% (light blue, inner) and 95% (dark blue, outer) 
confidence regions for the best-fit values of ao (in units of 
ms~ 2 ) and a for the MG11 and FS11 data sets (plotted on 
different sets of axes). It is apparent that the quality of the 
data makes it difficult to distinguish a = 4 (the prediction of 
MOND) from lower values; this difficulty is reflected in the 
wide range of slopes measured in other studies [7, 9, 10, 19- 
30]. 



MG11. We plot a Gaussian with the resulting mean 
(log[ao] = —9.85) and variance (a = 0.38 dex) as the blue 
dashed line in Fig. 1. The width of our Monte Carlo- 
generated distribution is much larger than that of the 
data values themselves (regardless of whether the latter 
width is taken to be 0.24 or 0.28 dex), so the data some- 
how show less scatter than would be expected from the 
errorbars. MGll's claim that the scatter in the data is 
explained by observational uncertainty alone overlooks 
this statistical feature. 

Performing the same analysis on the FS11 data, we 
find that a Gaussian fitted to the calculated a values has 
a = 0.33 dex. Using a Monte Carlo procedure as above 
to determine the most likely Gaussian that fits the data 
with the given errorbars, we find log(ao) = —9.82 and 
a = 0.37 dex. The FS11 data therefore show increased 
scatter as compared with the MG11 data, and although 
the scatter is still not completely accounted for by obser- 
vational uncertainty, there is less of a discrepancy than 
for the MG11 sample. 

Fitting the BTFR. — We determine the values of the 
BTFR parameters (ao and the slope, a) by minimizing 
the x 2 statistic that accounts for errors in both log(Mb) 
and V{, given in this case [47] by 



E 



[log(Mb) - alog(V{) + log(Ga )Y 



log(A/ h 



a 2 a? 

log(V f ) 



(3) 



where (Ji og (y t ) is determined by the asymmetric error bars 
on log(Vf ) that follow from assuming Gaussian uncertain- 
ties for Vf . This expression does not account for intrinsic 
scatter, since MOND predicts that it should be precisely 
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best-fit parameters 




ao free, a = 4 fixed 


both ao and a free 


data 


a x 10 10 (ms~ 2 ) 


x 2 


a x 10 10 (ms~ 2 ) 


log(a /[ms 2 ]) 


a 


x 2 


MG11 (47 galaxies) 
FS11 (58 galaxies) 


1.27 ±0.09 
1.33 ± 0.08 


65.0 
137.8 


0.761^ 
0.08t°;£ 


-10.1 ± 1.1 
-11.1 ± 0.8 


3.96 ±0.23 
3.75 ±0.17 


64.9 
133.0 



TABLE I: BTFR parameters that minimize the value of \ 2 [ E q. (3)] for the MG11 and FS11 data sets. The effect of the disk 
rotation factor described below Eq. (1) has been included in the calculation of ao. When a is left free, it is more natural to 
state the constraints on ao in terms of its logarithm, but we also give confidence intervals for ao itself for comparison with the 
fixed-slope case. The lowest \ 2 values we find are much larger than those found in MG11 (around 44). 



The parameter values that minimize Eq. (3) are given 
in Table I. The uncertainties on these values were deter- 
mined by the limits of the 68% confidence regions de- 
fined by x 2 < Xmin + Ax 2 j where A% 2 = 1-0 for a one- 
parameter fit or 2.3 for a two-parameter fit (see, e.g., 
Ref. [31]). Observe that the lowest x 2 value we can ob- 
tain, 64.9, is far from the —44 stated by MG11. 

Examining the confidence regions for the two- 
parameter fits (shown in Fig. 2), we find a broad range 
of slopes that can provide acceptable agreement with the 
data. A value of a — 4 is contained within the 95% con- 
fidence regions for both data sets, but these regions also 
cover slopes down to 3.6 for MG11 and 3.5 for FS11, a 
large range that reflects the variety of slopes measured 
in other studies (which we describe below). Either data 
set can only weakly discriminate between the MONDian 
slope (a = 4) and one bearing no connection to MOND 
(a ^4). 

Other measurements of the BTFR. — There are several 
studies that provide estimations of the BTFR slope from 
different galaxy samples, utilizing a wide variety of tech- 
niques for measuring Vf and Mb for each galaxy. Some 
find that slopes —4 best describe their data [3, 4, 7, 9, 19], 
while others find values roughly midway between 3 and 
4 [20-25] or lower still [10, 26-28], or alternately find a 
sensitive dependence on the methods used to convert lu- 
minosity to stellar mass [29, 30]. Useful summaries of 
these measurements are contained in Refs. [10] and [28]. 
Many of these studies stand in contrast with the results 
of MG11. Also, there are indications that if the mass of 
neutral hydrogen in a galaxy is used instead of baryonic 
mass, the resulting Tully-Fisher relation strongly favours 
Newtonian gravity over MOND [32]. 

Various arguments exist as to why some of the above 
measurements might be flawed. For example, Ref. [4] 
claims that the data used by Ref. [10] are not normally- 
distributed and hence cannot be used in a least-squares 
fit, and also that the rotation curves of Ref. [23] might 
not be extended enough to supply accurate values for Vf. 
Meanwhile, Ref. [10], conjectures that MGll's selection 
criteria of flat rotation curves could bias the determina- 
tion of a toward higher values, and Ref. [23] highlights 
issues surrounding the assembly of galaxy samples from 
several different observational surveys. 



When possible sample selection biases and differences 
in observational and data analysis techniques are taken 
into account, it becomes extremely difficult to rank the 
above studies in terms of accuracy of their estimations 
of the BTFR. It is clear, however, that the question of 
"the" slope of the relation is far from settled, certainly 
not at the level needed to test a theory predicting a value 
between 3 and 4. 

Cold dark matter and the BTFR.— Figure 2 of MG11 
plots that paper's galaxy data along with the MOND- 
predicted BTFR, Eq. (2), and a line taken to represent 
the BTFR anticipated in the CDM scenario. What is 
used for the latter is a cosmologically-motivated scaling 
relation (see, e.g., Ref. [33]), which has M h oc , un- 
der the assumptions that, for each galaxy, Mb is equal 
to the cosmic baryon fraction times the virial mass and 
V[ is equal to the rotation velocity at the virial radius. 
Neither of these assumptions is well-motivated, since the 
baryon fraction has been seen to vary from galaxy to 
galaxy (Ref. [23] provides a recent example of this), and, 
as mentioned by Ref. [10], the theory of Navarro- Frenk- 
White profiles leads us to believe that Vf is not repre- 
sentative of Kir- Ref. [4] takes steps towards improv- 
ing these assumptions by introducing extra fitting factors 
into the relations between the virial and observed quanti- 
ties, but still overlooks a fundamental problem with this 
approach. 

The essence of the problem is that galaxies are complex 
objects, with individual histories and properties deter- 
mined by feedback mechanisms, "gastrophysics," and a 
whole host of processes that are not yet fully understood. 
Therefore, a first-principles calculation of an analytical 
relationship between masses and rotation velocities of 
galaxies that is expected to hold in general is simply not 
possible in the context of CDM — a complete numerical 
calculation of the detailed physics that would affect such 
a relationship is still beyond our reach. Semianalytic 
treatments (e.g. [34]) can make some progress, but in 
reality, these efforts tend to be tuned to match observa- 
tions of the Tully-Fisher relation, rather than providing 
predictions. In this light, the fact that the supposed pre- 
diction of ACDM is nowhere near the data in MGll's 
figure 2 is misleading — there are certainly published es- 
timates for the BTFR, obtained through techniques like 



4 



halo abundance matching [35] , as well as models for feed- 
back processes [36, 37], that show tentative agreement 
with current data. Rcf. [4] describes how these studies 
are imperfect in various ways, but this is only to be ex- 
pected, as our understanding of galactic physics is not 
yet complete. 

Discussion. — MGll's conclusion, that the scatter- free 
BTFR describing a certain sample of galaxies strongly 
favours MOND over ACDM, is a bold one, and should 
be evaluated carefully and cautiously. We have raised 
concerns about the scatter of the data sample, the data's 
preference (or lack thereof) for a certain BTFR slope, the 
wide variety of measured slopes found in the literature, 
and what MG11 claims the theory of CDM has to say 
about the BTFR. These concerns are supplemented by 
those of Rcf. [38], which highlights how the measured 
baryonic masses in MG11 do not include the ionized gas 
content of galaxies; models that take this into account 
could potentially provide reasonable agreement both with 
the observed data and with recent N-body simulations of 
dark matter halos [39]. 

MOND, treated as a phenomenological description 
of galaxy-scale physics (particularly, galaxy rotation 
curves [40, 41] and the properties of dwarf galaxies [42]), 
is apparently quite successful. However, the tenets of 
MOND are challenged in the same regime by models of 
the globular cluster NGC 2419 [43, 44] and fits to the 
"neutral hydrogen Tully-Fisher relation" [32]. As well, 
the scalings and value of ao are unsurprising when put in 
the context of familiar facts about galaxies and cosmol- 
ogy [45]. 

Due to the ambiguity present in current data and in our 
picture of the behaviour of galaxies, the BTFR provides 
no new advantage for MOND over CDM. Nevertheless, as 
the data improve and systematics become less significant, 
the BTFR could begin to provide useful information. In 
particular, more precise measurements of a could serve 
as a guideline with which new models for supernova feed- 
back would have to agree, and these measurements could 
therefore assist in the development of a full description of 
the baryonic processes that dominate deep within CDM 
haloes. It is also worth exploring variants of the BTFR, 
using other mass tracers like neutral hydrogen [32] or 
considering the relation at different radii [46] (although 
it is nontrivial to define a suitable radius consistently for 
gas-dominated galaxies) . If supplemented with other ob- 
servational probes, the BTFR could well provide a unique 
window into the realm of galactic physics. 
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